The Reconstruction of Pt(lll) and 
Domain Patterns on Close-packed Metal Surfaces 
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We have studied the reconstruction of Pt(lll) theoretically using a two-dimensional Frenkel- 
Kontorova model for which all parameters have been obtained from ab initio calculations. We 
find that the unreconstructed surface lies right at the stability boundary, and thus it is relatively 
easy to induce the surface to reconstruct into a pattern of FCC and HCP domains, as has been 
shown experimentally. The top layer is very slightly rotated relative to the substrate, resulting in 
the formation of "rotors" at intersections of domain walls. The size and shape of domains is very 
sensitive to the density in the top layer, the chemical potential, and the angle of rotation, with a 
smooth and continuous transition from the honeycomb pattern to a Moire pattern, via interlocking 
triangles and bright stars. Our results show clearly that the domain patterns found on several 
close-packed metal surfaces are related and topologically equivalent. 

PACS numbers: 68.35.Bs, 61.72.bb, 68.35.Md, 68.55.-a 
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Due to its enormous importance as a catalyst, Pt(lll) 
is one of the most widely studied surfaces. Under "nor- 
mal" conditions, it has the flat topography expected of a 
bulk-truncated face-centered-cubic (FCC) (111) surface. 
However, experiments have shown that one can induce 
the surface to reconstruct into either a honeycomb struc- 
ture or a pattern of interlocking triangles - for example, 
by heating the surface above 1330 K J5J], or by placing 
it in a supersaturated Pt vapor J2|. The reconstructed 
structure is comprised of domains where the bulk FCC 
stacking sequence is retained, alternating with domains 
where the surface atoms instead occupy hexagonal-close- 
packed (HCP) sites, sitting directly above atoms two lay- 
ers below. 

This is one of a family of similar reconstructions, 
formed by a tessellation of FCC and HCP domains, seen 
on Au(lll) as well as various heteroepitaxial systems on 
the (111) faces of FCC metals, and the structurally sim- 
ilar (0001) faces of HCP metals §, |, §, §]. These struc- 
tures have attracted a great deal of interest, especially 
because they can be used as templates for growing or- 
dered arrays of nanoparticles J?], |[ Possible appli- 
cations for such nanostructures include nanoelectronics, 
information storage, and nanoscale chemical reactors. To 
design such nanostructures, it is desirable to understand 
the factors controlling the geometry and spacing of the 
reconstruction patterns. 

In this paper, we study the structure of the Pt(lll) 
surface theoretically. We show that the unreconstructed 
Pt(lll) surface is in fact teetering right at the brink of 
a domain of stability. Thus, slight changes in the envi- 
ronment can trigger a reconstruction, whose periodicity 
and geometry are very sensitive to various extrinsic and 
intrinsic parameters. In addition to obtaining excellent 
agreement with the structures reported experimentally 
for Pt(lll), we also observe most of the structures re- 



ported experimentally for other systems. 

The presence or absence of such reconstructions in- 
volves a very delicate balance between various contribu- 
tions to the total energy, and it is therefore desirable 
to have as accurate a description of interatomic interac- 
tions as possible. Unfortunately, the very large unit cells 
of the reconstructed surfaces make it unfeasible to per- 
form a fully ab initio calculation. However, one can use 
ab initio calculations to parametrize a model; this is the 
approach we will follow. 

The driving force for the reconstruction of these close- 
packed surfaces is tensile surface stress, i.e., surface 
atoms would like to be closer together than the bulk 
nearest- neighbor (NN) spacing a [^5|. This tendency is 
however opposed by the fact that (in general) it costs en- 
ergy when: (i) surface atoms lose registry with the sub- 
strate (ii) extra atoms must be provided to increase the 
surface density. These factors are incorporated in the 
Frenkel-Kontorova (FK) model, which has been widely 
used to study commensurate- incommensurate transitions 
]l0| . In its original form, the model consists of a 
one-dimensional chain of atoms connected by harmonic 
springs of equilibrium length b a, sitting in a sinusoidal 
substrate potential. This has the advantage of being ex- 
actly solvable, and can serve as a useful guideline for 
studying structural stability. 

However, a proper description requires a generalized 
version of the FK model, with a two-dimensional (2D) 
layer of atoms interacting via a more realistic anharmonic 
potential V ss , and sitting in a 2D potential due to bulk 
atoms, V s b, with competing minima at FCC and HCP 
sites. The Hamiltonian is given by [ fll) : 
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where i runs over all NN bonds of length U between sur- 
face atoms, and j runs over all atoms at positions in 
the surface layer. TV is the total number of atoms, and T 
is a chemical potential that contains information about 
the energy required to incorporate an atom into the sur- 
face layer; this depends on where the atom comes from 
(bulk, step edge, adatom, etc.) 

To obtain V ss , V s b and V, we have performed ab initio 
density functional theory calculations to study the ener- 
getics of compressing the surface layer, making surface 
stacking faults, and extracting atoms from various sites. 
We use the PWSCF package Q , with a plane wave basis 
with a cut-off of 20 Ry, ultrasoft pseudopotentials [^3| , 
and the local density approximation to the exchange- 
correlation potential. Surface calculations have been car- 
ried out using supercells with 9 layers of atoms, separated 
by a vacuum of 6-layer thickness. A Monkhorst-Pack grid 
corresponding to 27 k-points in the irreducible part of the 
surface Brillouin zone is used to sample reciprocal space 
for calculations using a (1 x 1) surface cell; the grid is var- 
ied commensurately when using larger surface cells. We 
obtain the bulk lattice constant as arj = 3.92A, and find 
that the first interlayer spacing 0I12 is slightly expanded 
by 0.3% relative to the bulk interlayer spacing, and the 
next two interlayer spacings 0I23 and ^34 are slightly con- 
tracted by 0.55% and 0.18% respectively. We obtain the 
(unreconstructed) surface energy 7 and surface stress a 
as 9.13 and 29.45 mRy/A 2 respectively. These results 
are in good agreement with experiments and previous 
calculations [l5|, [r| . 

The trickiest part of the parametrization concerns V ss , 
the interaction between surface atoms. Earlier authors 
(who studied other systems), either used "physically rea- 
sonable" parameters |T^] , or did calculations to study the 
energetics of a monolayer of atoms on a jellium (T^|; the 
ambiguity in the latter approach arises from the uncer- 
tainty in the choice of jellium density. We have chosen 
instead to calculate the variation in the surface stress 
when compressing a slab of atoms |^6| . The distance d\2 
is allowed to relax, and it is therefore a good approxima- 
tion to assume that the variation in surface stress comes 
entirely from the surface-surface bonds that we are in- 
terested in. We use these results to parametrize a Morse 
potential V ss = A${1 — exp[— A\(a — b)]} 2 ; we obtain Aq 
= 60.2 mRy, A x = 2.062 A" 1 , and b = 2.638 A, i.e., 
surface bonds would like to shorten their length by 4.7% 
from the bulk NN distance a = a /^2 = 2.77 A. 

To obtain V s b, we study surface stacking faults. We 
find that the energy cost (per surface atom) of having the 
surface layer occupy HCP, top and bridge sites instead of 
the most favored FCC site is 5.05, 11.85 and 5.89 mRy 
respectively. These values are then used to expand V s b in 
a 2D Fourier series, using the first two shells of reciprocal 
lattice vectors for the 2D surface lattice |l8|, Fig. 1 
shows V s b for a line cutting through the surface along 
the [112] direction. Note that though the HCP site is a 
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FIG. 1: The inset shows a top view of the Pt(lll) surface, 
while the main graph shows the surface-bulk potential V s b 
along the line PQRS marked in the inset. 

local minimum, it lies considerably above the FCC site 
in energy. 

Finally, to obtain T, we use our ab initio results for 
the bulk cohesive energy, the surface energy, the adatom 
adsorption energy and the energy needed to detach an 
atom from a kink site at a step edge |^0) to obtain Ti, = 
50.56 mRy for bulk atoms, T a = -67.12 mRy for adatoms, 
and Tfc = 50.4 mRy for kink atoms. 

To see how stable the unreconstructed surface is, one 
can map the 2D problem onto a ID FK model |ll[, and 
evaluate the dimensionless parameter R = (37ra)(|(7 — 
7)/8v / 2fcVF, where k is the spring constant for NN sur- 
face bonds and W is the amplitude of V s b along the 
zigzag line connecting adjacent FCC and HCP sites. If 
R > / < 1, the surface will / will not reconstruct. Mans- 
field and Needs Jll[] have obtained R = 0.73 for Pt(lll) 
and R = 0.44 for Au(lll), i.e., they predict that unre- 
constructed Au(lll) should be more stable than Pt(lll), 
whereas it is well known that Au(lll) reconstructs even 
at low temperatures. We believe that the deficiency lies 
not with their model but with the values they used to 
obtain R. Upon approximating V ss by a harmonic poten- 
tial in the region of interest, we obtain k—370 mRy/ A 2 . 
Using W = 5.89 mRy, we get R = 0.998, i.e., we find 
that the unreconstructed surface is only just stable. This 
explains why it is relatively easy to induce the surface 
to reconstruct by increasing the temperature (reducing 
W) or by placing the surface in a supersaturated vapor 
(effectively increasing the numerator of R). 

For a more precise determination of structural stabil- 
ity, we work with the full 2D Hamiltonian of Eq. 1. Wc 
study the variation in the surface energy 7 as a function 
of Ap, the excess density in the surface layer (relative to 
the unreconstructed case), and 9, the angle between the 
unit vectors of the top layer and the substrate. (Though 
optimizing 9 makes only a very small difference to 7, we 
will show below that allowing for 9 7^ is crucial to ex- 
plain the "rotors" at intersections of domain walls.) In 
all cases, we start with a uniformly compressed layer of 
surface atoms, and obtain the atomic positions that min- 
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FIG. 2: Results from the 2D Frenkel-Kontorova model for 
isotropic reconstruction, (a) and (b) show how A7, the 
change in the surface energy (relative to the unreconstructed 
surface) and 6, the angle of rotation of the top layer, depend 
on the excess density Ap; (c) and (d) show how the opti- 
mal excess density Ap m i n and the periodicity I vary with the 
chemical potential V. 



imize 7, using a conjugate gradient algorithm. 

For Pt(lll), we find that isotropically compressed 
structures, with threefold symmetry and sixfold atomic 
coordination everywhere, are always lower in energy than 
the uniaxial "stripe" patterns or structures with point 
dislocations. Fig. 2(a) shows how A7, the difference 
in 7 for the reconstructed and unreconstructed surfaces, 
varies with Ap. Each point is also optimized with respect 
to 8; Fig. 2(b) shows that 8 varies approximately linearly 
with Ap and is very small. If the surface reconstructs, 
A7 will have a minimum at a non-zero value Ap m i n . 
Under normal conditions, extra atoms are obtained from 
the bulk or step edges, and T ~ T^. It is clear from the 
figure that under these conditions, the surface will not 
reconstruct, which agrees with experiment. However, if 
adatoms are available for incorporation in the surface 
layer, T is lowered, up to a minimum of r a , and the sur- 
face then reconstructs. Fig. 2(c) shows how Ap min varies 
as T varies between r o and Tf,. The surface reconstructs 
only if r < 36 mRy. As T decreases and Ap increases, 
the periodicity / of the reconstruction decreases; Fig 2(d) 
shows how I varies with T. 

Our most striking results are obtained upon examining 
the domain patterns using a simple technique to obtain 
the surface corrugation and thus simulate STM images. 
The height of each atom in the unit cell is obtained using 
a 2D Fourier expansion similar to that used for V s b, ex- 
panding about the heights at FCC, HCP, bridge and top 
sites. From our ab initio calculations on surface stack- 
ing faults, we obtain these heights as 0, 0.03, 0.29 and 
0.06 A respectively. These are smaller than the corru- 
gations measured experimentally However, we find 




FIG. 3: Simulated STM images of the Pt(lll) surface, as 
Ap is varied from 2.9% to 21%. Individual atoms are shaded 
according to their height; lighter atoms are higher. In (a) 
the top layer is aligned with the substrate, in (b) to (f), the 
angle betwen the top layer and substrate has been optimized, 
(a) and (b) are the honeycomb, (c) and (d) are threefold 
whorls (e) is the bright star, and (f) is the Moire pattern. 
The black/ white line in each image has a length of 50 A 



that ab initio simulations of the reconstructed structure 
(for small cell sizes accessible to computation) support 
our smaller values ^0|. As has been reported recently 
for Au(lll), it appears that the STM exaggerates the 
value of the corrugation pl[ ; in any event, we are pri- 
marily interested in the domain patterns, rather than in 
the absolute heights. 

A sequence of simulated STM images is shown in 
Fig. 3. In all six images, black, dark grey, light grey 
and white areas correspond to regions where atoms sit 
at FCC, HCP, bridge and top sites respectively. A non- 
linear grey scale has been used to show the domain walls 
(bridge sites) clearly. Fig. 3(a) shows the pattern ob- 
tained when Ap=2.9%, and 8 = 0. Large hexagonal 
FCC domains are separated from narrow HCP domains 
by domain walls that form a honeycomb network. Alter- 
nate vertices of the hexagons are extremely bright; atoms 
here sit at top sites. Between these bright vertices, there 
is a "three-pointed star" similar to that obtained in an 
earlier simulation of the nucleation of the reconstruction 
p2|. Fig. 3(b) shows the structure obtained for the same 
Ap, but with 8 at its optimal value of 0.14°. It is evident 
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from a comparison of the two images that breaking the 
symmetry by even this tiny amount has a big impact on 
the structure. The bright vertices are now transformed 
into "bright rotors" , which are indeed observed experi- 
mentally j2j|. The similarity between Fig 3(b) and the 
experimental STM images |2j| is striking. Our find- 
ing that the rotors arise from a rotation of the top layer 
relative to the substrate explains why, in a given region, 
experiments tend to find either all clockwise or all an- 
ticlockwise rotors. The honeycomb pattern occurs in a 
region where I is very sensitive to T; we believe this ex- 
plains why STM images usually show irregular hexagons. 

As r is lowered, the surface layer densities further. 
Figs. 3 (c) and (d) show the structures obtained when 
Ap is increased to 4.0% and 6.8% respectively. The hon- 
eycomb transforms into a pattern of interlocking wavy 
triangles. Such structures have been observed on dense 
terraces of Pt(lll) pj|. This structure has also been seen 
for Na/Au(lll) § and for 3 ML of Cu/Ru(0001) @. As 
r is lowered and Ap increased still further, the wavy do- 
main walls straighten out, resulting in the "bright star" 
pattern shown in Fig. 3(e), for a value of Ap— 10.2%. 
This pattern has not been seen on Pt(lll) as it requires 
a very low T (high adatom coverage) ; however it is seen on 
other surfaces, e.g., Ni/Ru(0001) §. Ap min for Pt(lll) 
has a maximum value of 12% (when T = T a ); it is 
nevertheless instructive to see what happens for larger 
Ap. The triangular domains of the bright star transform 
into hexagons, resulting in the Moire pattern shown in 
Fig 3(f). The density distribution on the surface is now 
almost uniform. This pattern is also observed experimen- 
tally, e.g., for 4 ML of Cu/Ru(0001) §. 

The sequence of images shown in Fig. 3 makes it clear 
that all these structures "morph" smoothly into one an- 
other and are topologically equivalent. Progressively 
larger values of Ap can be favored either by lowering T [as 
for Pt(lll)], or by increasing the surface stress - e.g., by 
depositing more overlayers in a heteroepitaxial system. 
Indeed, a similar progression of structures with overlayer 
thickness has been observed for the Cu/Ru(0001) system 
H Our results show that the same physics is oper- 
ating in all the systems mentioned above, and it is only 
slight changes in parameters that are responsible for the 
various domain patterns observed on different surfaces. 

To summarize: we have shown that the unrecon- 
structed Pt(lll) surface is only just stable, and can thus 
be easily induced to reconstruct. We have confirmed that 
Pt(lll) will reconstruct if the chemical potential is low- 
ered by the presence of a large number of adatoms. In 
addition to the methods tried to date, we suggest that 
it should also be possible to induce the surface to recon- 
struct in an electrochemical environment or by depositing 
alkali metals on it. The domain patterns we obtain for 
Pt(lll) arc in excellent agreement with experiment. We 
have also shown that slight variations in conditions can 
lead to many of the other domain patterns seen on sim- 



ilar surfaces. We have shown that the periodicity of the 
reconstruction can be controlled by varying the chemical 
potential; this is important if one would like to use such 
surfaces as templates for growing ordered nanostructures. 

We thank S. de Gironcoli for providing the Pt ultrasoft 
pscudopotential, and X. Gonze, I.K. Robinson, H. Brunc 
and U.V. Waghmare for helpful discussions. 
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